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A  GEOACOUSTIC  MODEL  FOR  FINE-GRAINED,  UNCONSOLIDATED, 
CALCAREOUS  SEDIMENTS  (ARSRP  NATURAL  LABORATORY) 


OBJECTIVE: 

The  long-range  scientific  goal  of  the  Acoustic  Reverberation 
Special  Research  Program  (ARSRP)  is  to  understand  and  predict  ocean 
broadband  reverberation  at  low  frequencies  (50-1000  Hz)  over  long 
propagation  paths  (up  to  1000  km).  To  implement  this  goal,  the  concept 
of  Natural  Laboratories  in  the  Atlantic  and  Pacific  Oceans  (Tucholke  et 
al.,  1991)  was  developed  as  an  effective  means  for  studying  geological 
and  geophysical  variability  of  seabed  properties  that  determine  acoustic 
bottom/subbottom  reverberation.  This  memorandum  report,  in  support 
of  ARSRP  reverberation  modeling  in  the  Atlantic  area,  represents  a 
compilation  of  sediment  (geoacoustic)  properties  deemed  important  as 
input  parameters  to  acoustic  field  models. 

APPROACH: 

Because  the  volume  of  published  sediment  geoacoustic  data  for 
the  Atlantic  Natural  Laboratory  region  is  small,  the  database  was 
expanded  to  include  geoacoustic  data  from  other  ocean  basin  areas  in 
order  to  obtain  better  statistical  representations  for  each  parameter. 
The  Natural  Laboratory  is  located  east  of  Bermuda  on  the  flanks  of  the 
Mid-Atlantic  Ridge  and  has  bottom  depths  that  are  generally  shoaler 
than  the  calcium  compensation  depth.  It  is  assumed,  therefore,  that  the 
sediments  will  be  pelagic,  dominantly  fine-grained,  calcareous,  and  will 
have  generally  the  same  physical  characteristics  as  calcareous 
sediments  in  other  open  ocean  areas  (i.e.,  the  deposition  of  a  calcareous 
clay/ooze  is  essentially  the  same  for  all  deep-ocean  areas). 

The  geoacoustic  data  are  derived  primarily  from  extensive 
laboratory  and  shipboard  measurements  made  on  Deep  Sea  Drilling 
Project  and  Ocean  Drilling  Program  (DSDP/ODP)  cores  and  published  in 
the  DSDP/ODP  cruise  reports  (Table  1).  Compositionally,  the  sediments 
in  the  database  range  in  carbonate  content  from  30%  to  nearly  100^ 
(Fig.  1).  Lithologically,  they  are  primarily  oozes  and  chalks.  A 
substantial  amount  of  geoacoustic  data  in  the  DSDP/ODP  reports  is  not 
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Table  1.  DSDP  and  ODP  Sources  of  Sediment  Data 


DSDP  Site 

Leg 

Samples 

Location 

400A 

48 

30 

North  Atlantic 

401 

M 

22 

if 

ti 

516F 

72 

18 

South  Atlantic 

ODP  Site 

Leg 

SampJ&s 

Location 

628A 

101 

7 

Bahama  Bank 

630 

•t 

10 

ii 

ii 

632A 

" 

6 

Exuma  Sound 

633 

ii 

4 

it 

647 A 

105 

15 

Labrador 

Sea 

657  A 

108 

41 

Equatorial 

Atlantic 

657B 

t« 

5 

ti 

" 

662A 

ii 

12 

it 

ii 

663A 

ii 

8 

" 

" 

664B 

ii 

12 

ii 

it 

705A 

115 

6 

Mascarene 

Plateau 

706A 

11 

13 

ii 

ii 

706B 

ti 

12 

it 

ii 

707A 

ii 

16 

ti 

707B 

" 

14 

" 

" 

707C 

tt 

6 

it 

it 

708A 

II 

28 

H 

it 

709A 

" 

43 

it 

n 

709B 

" 

31 

it 

tt 

710A 

II 

2 

it 

" 

710B 

" 

2 

ii 

it 

71 1 A 

" 

11 

it 

" 

711B 

II 

3 

ii 

ii 

712A 

It 

9 

tt 

it 

713A 

" 

1 

it 

" 

714A 

" 

41 

" 

" 

714B 

II 

18 

it 

" 

728A 

117 

10 

Arabian 

Sea 

730A 

" 

12 

ii 

it 

731 A 

ii 

21 

ii 

ii 

744A 

119 

12 

Kerguelen  Plateau 

758A 

121 

47 

Ninetyeast  Ridge 

760A 

122 

16 

Wombat 

Plateau 

761 B 

•i 

35 

if 

it 

762B 

ti 

24 

Exmouth 

Plateau 

762C 

♦i 

40 

M 

it 

763A 

ii 

23 

It 

it 

764A 

ii 

7 

Wombat 

Plateau 

764B 

ii 

2 

ii 

” 

765  B 

123 

5 

Argo  Abyssal  Plain 

766A 

ii 

8 

Exmouth 

Plateau 

768  B 

124 

12 

Sulu  Sea 

770B 

it 

4 

Celebes  Sea 

787B 

126 

5 

Izu-Bonin 

Basin 

790C 

” 

5 

Minami-Sumisu  Basin 

791B 

•i 

16 

ii 

it  it 

793A 

ii 

9 

Izu-Bonin 

Basin 

793  B 

ti 

11 

tt  it 

it 

2 


included  in  the  database  for  the  following  reasons:  (1)  geoacoustic 
measurements  with  no  accompanying  carbonate  content  information 
were  excluded.  The  importance  of  clearly  defining  sediment  type  (i.e., 
ooze  vs.  clay)  in  relationship  to  the  various  properties  is  important  for 
the  extrapolation  of  these  properties  to  unsampled  or  poorly  sampled 
areas  of  the  Natural  Laboratory;  and  (2)  graphical  presentations  of 
analytical  results  were  usually  not  detailed  enough  to  be  usable. 

Results  from  extensive  studies  on  the  acoustic  and  physical 
properties  of  deep-sea  carbonate  sediments  reported  by  Johnson  et  al. 
(1977),  Milholland  et  al.  (1980),  Mayer  (1979),  and  Hamilton  et  al. 
(1982)  are  not  included  in  this  compilation  but  are  summarized  along 
with  this  study  in  Table  2.  The  study  by  Hamilton  et  al.  (1982)  consists 
of  323  samples  and  includes  the  work  of  Johnson  et  al.  (1977)  and 
Mayer  (1979).  The  samples,  which  are  from  the  Ontong-Java  Plateau 
and  the  eastern  equatorial  Pacific,  represent  the  upper  9  m  of  sediment; 
roughly  half  of  the  samples  occur  in  the  top  45  cm  of  sediment. 
Milholland  (1978)  made  a  systematic  study  of  269  DSDP  core  samples 
representing  a  complete  ooze-chalk-limestone  sequence  (1262  m)  on 
the  Ontong-Java  Plateau  (sites  288  and  289). 

Relative  to  the  geoacoustic  properties  previously  mentioned, 
compressional  wave  attenuation,  shear  wave  attenuation,  and  shear 
wave  velocity  are  sparsely  represented  in  the  literature.  Consequently, 
in  order  to  obtain  reasonable  estimates  of  these  properties,  all  available 
data  were  used  that  related  to  fine-grained,  or  predominantly  fine¬ 
grained,  sediments.  In  cases  where  the  sediment  is  designated  as 
turbidite,  sand-silt-clay,  or  mud,  it  is  assumed  that  there  is  a  substantial 
fine-grained  component.  In  abyssal  plain  areas,  for  example,  turbidite 
layers  are  usually  thin  relative  to  the  clay-ooze  intervals  that  separate 
them.  In  the  case  of  shear  wave  attenuation,  so  little  data  are  available 
that  some  measurements  made  in  sandy  sediments  are  included  in  the 
database. 

The  data  in  this  report  are  presented  as  (1)  tables  that  summarize 
various  statistical  parameters  for  the  sediment  geoacoustic  properties 
previously  mentioned,  (2)  cross  plots  that  illustrate  the  relationships 
between  various  properties,  and  (3)  a  geoacoustic  model  based  on  1  and 
2.  Curve  fits  to  the  data  are  provided  by  regression  equations  (data 
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permitting).  Regression  analysis  represents  a  means  of  quantifying  the 
relationships  among  sediment  properties.  The  equations  can  be  used  to 
predict  values  for  a  desired  property  when  values  of  another  physically 
related  property  are  known. 

When  summarizing  analytical  measurements  (mean  values, 
regression  equations,  etc.)  from  a  variety  of  sources  (Table  1),  errors 
are  inherent  due  to  differences  in  measurement  technique,  instrument 
calibration  and  precision,  sample  disturbance,  sampling  bias,  etc. 
Nevertheless,  Table  2  shows  that  the  results  of  this  study  compare  well 
with  similar  studies.  The  difficulty  in  measuring  shear  waves  and 
determining  attenuation  in  sediments,  including  the  fact  that  the 
mechanisms  governing  these  parameters  are  poorly  understood,  make 
them  the  least  reliable  of  all  the  geoacoustic  parameters. 

DISCUSSION  OF  DATA: 

General: 

Measured  geoacoustic  values  extracted  from  the  DSDP/ODP 
literature  are  summarized  in  Table  3.  The  wide  range  in  values  is  due 
largely  to  the  factors  mentioned  above  but  also  because  the  individual 
values  represent  a  wide  sampling  range  within  the  sediment  column. 
Chalks  and  oozes  are  listed  separately  because  they  are  (usually) 
lithologically  distinct  and,  therefore,  have  different  ranges  of  values  for 
the  same  properties.  DSDP  lithologic  nomenclature  conventions 
(Winterer  et  al.,  1973)  define  oozes  as  soft  and  readily  deformed  using 
the  finger  or  broad  spatula  blade.  Chalks,  which  are  also  soft,  are  less 
easily  deformed  (i.e.,  chalks  are  partly  indurated  oozes). 

Figures  2,  3,  and  4  show  reasonably  good  empirical  relationships 
between  compressional  wave  velocity,  bulk  density,  and  porosity. 
Within  the  limits  defined  in  Table  3,  compressional  velocity,  shear 
velocity,  and  bulk  density  all  increase  in  magnitude  with  depth  below 
the  seafloor,  while  porosity  decreases  (Figs.  5,  6,  7,  and  8).  Curve  fits 
for  these  depth  plots  appeared  reasonable  only  in  the  case  of  bulk 
density  (Fig.  6)  although  a  fit  is  provided  for  shear  wave  velocity  as 
well  (Fig.  8).  No  relationship  was  established  between  shear  wave 
velocity  and  the  other  properties. 
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Compressional  velocity  (Fig.  9),  density  (Fig.  10),  and  porosity  (Fig. 
11)  appear  to  be  independent  of  carbonate  content.  Hamilton  et  al. 
(1982)  have  noted  that  as  carbonate  content  exceeds  75%  velocities  may 
increase  due  to  the  presence  of  coarse  sand-sized  particles  (e.g., 
foraminifera  tests).  This  database  consists  of  deep-water  oozes  (2000- 
4000  m)  containing  little  or  no  coarse  material.  Consequently,  no 
increase  in  velocity  is  observed  at  carbonate  percentages  above  75%. 

Hamilton  et  al.  (1982)  also  noted  a  correlation  between  bulk 
density  and  percent  carbonate  in  sample  suites  where  changes  in 
carbonate  content  are  balanced  by  changes  in  biogenic  silica  content. 
If,  on  the  other  hand,  the  changes  in  carbonate  content  are  balanced  by 
changes  in  pelagic  clay  content  (similar  grain  density  as  calcium 
carbonate)  no  usable  relationship  exists.  Samples  containing  biogenic 
silica  were  omitted  from  the  database  forming  this  report.  Therefore,  it 
follows  that  no  correlation  was  observed  between  bulk  density  and 
percent  carbonate. 

Compressional _ Wave _ Meloeity* 

Urmos  and  Wilkens  (1993)  observed  that  the  presence  of  chalks 
coincided  with  greater  scatter  in  their  data.  Likewise,  the  scatter  in 
Figures  2-8,  which  includes  both  ooze  and  chalk  data,  is  attributed 
mainly  to  chalks  as  demonstrated  when  these  lithologies  are  separated. 
In  Figure  12,  for  example,  most  of  the  oozes  cluster  above  250  m, 
whereas  the  chalks  (Fig.  13)  are  widely  scattered  throughout  the 
sediment  column. 

The  observation  that  the  oozes  occur  predominantly  within  the 
top  250  m  of  sediment  suggests  that  soft,  pelagic,  carbonate  oozes 
transition  to  chalks  at  roughly  this  depth.  Urmos  and  Wilkens  (1993) 
place  the  ooze-chalk  transition  nearer  to  300  m.  These  observations  are 
consistent  with  Schlanger  and  Douglas  (1974)  who  defined  a  shallow- 
burial  realm  where  gravitational  compaction  is  the  dominant  process 
and  cementation  is  a  subordinate.  The  fact  that  chalks  are  found 
scattered  throughout  the  entire  depth  range  in  Figure  13  indicates  that 
the  ooze-chalk  transition  is  not  strictly  time  and  depth-of-burial 
dependent.  Cementation  can  be  enhanced  in  the  shallow-burial  realm 
for  a  variety  of  reasons,  resulting  in  the  development  of  chalk  above 
250  m.  Rapid  and  extensive  diagenesis  is  characteristic  of  carbonate 
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sediments  of  the  Bahama  Banks  area.  Exuma  sound  samples  from  ODP 
Site  633  account  for  all  but  one  of  the  chalks  that  occur  above  50  m  in 
Figure  13. 

The  scatter  exhibited  in  Figure  13  is  an  indication  of  the  variation 
in  lithification  and/or  cementation  associated  with  chalks  which,  in  turn, 
may  cause  uncertainty  in  their  identification.  A  number  of  the  "oozes" 
in  Figure  12,  for  example,  plot  well  outside  the  general  cluster  of  points. 
These  outlying  velocities  are  more  appropriate  for  chalks.  Similarly, 
some  "chalks"  (Fig.  13)  have  velocities  that  are  indicative  of  oozes.  A 
simple  linear  fit  can  be  applied  to  the  ooze  data  above  250  m  (Fig.  14). 
The  fit  is  slightly  improved  in  Figure  15  by  assuming  that  velocities  of 
1.60  km/sec  or  above  are  indicative  of  chalk  and  that  anything  less  is  an 
ooze.  Thus,  some  of  the  data  in  Figure  14  are  excluded  from  Figure  15 
while  some  of  the  data  in  Figure  13  are  included.  In  Figures  16  and  17 
bulk  density  and  porosity  have  been  "filtered"  in  a  similar  manner. 

The  regression  equations  in  Figures  15-17  are  useful  for  modeling 
the  upper  250-300  m  of  pelagic  sediment  if  one  assumes  that  there  are 
few  if  any  chalk  layers  within  this  interval.  This  assumption  may  not 
be  unrealistic.  Sonic  velocity  logs  recorded  under  laboratory  conditions 
(see  various  DSDP/ODP  reports),  in  addition  to  laboratory  measurements 
made  on  discrete  samples,  typically  represent  the  upper  few  hundred 
meters  of  sediment  as  having  relatively  uniform  compressional  wave 
velocity  profiles  (i.e.,  no  large  velocity  spikes). 

Urmos  and  Wilkens  (1993)  have  reported  a  discrepancy  between 
sediment  properties  measurements  made  in  the  laboratory  and  in  situ 
values.  They  attribute  this  to  the  hydraulic  rebound  of  pore  water 
when  sediment  is  transferred  from  ambient  conditions  to  laboratory 
conditions.  It  is  apparent  from  their  work  that  laboratory 

measurements  do  not  preserve  the  nature  of  the  velocity  gradients 
found  in  undisturbed  deposits.  As  a  result,  they  have  derived  an 

empirical  velocity  correction  factor  (AVp)  that  can  be  applied  to 
unconsolidated  calcareous  sediments  with  carbonate  contents  greater 
than  60%.  A  corrected  velocity  profile,  obtained  by  using  their 
technique,  is  shown  in  Figure  15. 
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Shear _ WaY.fi _ Velocity: 

Propagating  and  receiving  shear  waves  in  unconsolidated 
sediments  is  dnficult,  even  under  controlled  conditions.  For  this  reason 
and  those  given  in  the  earlier  discussion  of  approach,  one  could  question 
the  reliability  of  the  shear  wave  velocity  data  extracted  from  the 
DSDP/ODP  reports.  Although  the  scatter  in  the  data  is  substantial  in 
Figure  8,  a  linear  curve  fit  of  the  data  is  provided,  and  the  equation 
shown  is  used  to  compute  curve  6  in  Figure  18.  For  comparison,  eight 
other  curves  (some  commonly  referenced)  are  shown,  and  their 
associated  equations  are  listed  below.  The  equations  are  indexed 
according  to  the  number  given  to  each  curve  in  Figure  18. 

1.  Vs  =  78.89d0-3l2 

2.  Vs  =  116  +  4.65d  (where  d  =  0-36  m) 

Vs  =  237  +  1.28d  (where  d  =  36-120  m) 

Vs  =  322  +  0.58d  (where  d  =  120-650  m) 

3.  Vs  =  42.7064d0-3069 

4.  Vs  =  62.5223d0-3537 

5.  Vs  =  96.7977d°-29326 

6.  Vs  =  77.463  +  0.24079d 

7.  Vs  =  33.48  +  1.44d  +(-3.79xl0*3)d2  +  6.67xl0*6d3 
This  equation,  which  will  reproduce  the  curve 
shown  in  Figure  18,  is  based  on  the  relationships 
of  Bryan  and  Stoll  (1988)  that  require  void  ratio 
as  an  input.  In  the  case  of  curve  7,  void  ratios 
were  estimated  from  porosity  values  given  in 
the  geoacoustic  model  section  of  this  report. 

8.  Vs  =  75.4  +  1 .5  Id  +  (-5.47xl0*3)d2  +  8.35xl0'6d3 

9.  Vs  =  23.0d°-3.  The  initial  value  of  23.0  m/sec  is 
the  average  shear  wave  velocity  for  the  fine¬ 
grained  sediments  found  in  Table  1  from 
Richardson  et  al.  (1987).  The  power  of  0.3  with 
depth  (d)  is  from  Richardson  et  al.  (1991). 
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It  is  evident  in  Figure  8  that  most  of  the  shear  wave  velocities  in 
the  top  100  m  of  sediment  fall  within  25-125  m/sec  and  that  over  the 
full  depth  range  (300  m)  most  values  are  less  than  200  m/sec.  These 
values  are  significantly  lower  than  predicted  by  curves  1,  2,  4,  and  5  in 
Figure  18,  but  are  in  the  range  predicted  by  the  other  curves.  Curves  6 
and  9,  in  particular,  are  comparable  to  laboratory  measurements  made 
on  carbonate  sediments  by  Schultheiss  (1985)  and  Lavoie  and  Anderson 
(1991). 

Most  of  the  equations  above  predict  higher  surface  shear  wave 
velocities  than  are  usually  measured  in  muddy  sediments.  Direct  probe 
measurements  by  Richardson  et  al.  (1991),  for  example,  recorded 
velocities  of  20-25  m/sec  at  10-30  cm  and  40-50  m/sec  at  2  m  depth  as 
well  as  gradients  of  4-17  s1.  All  the  curves,  with  the  exception  of  6,  also 
predict  a  generally  rapid  increase  in  shear  wave  velocity  with  depth  for 
the  upper  50  m  of  sediment  followed  by  a  more  gradual  increase.  The 
laboratory  measurements  by  Schultheiss  (1985),  and  Lavoie  and 
Anderson  (1991)  do  not  show  any  significant  depth  dependence 
(maximum  sample  depth  =  270  m)  and,  in  this  respect,  their  data  agree 
more  closely  with  curve  6  in  Figure  18.  The  absence  of  any  significant 
gradient  in  their  data  indicates,  once  again,  that  velocity  measurements 
made  in  the  laboratory  do  not  reflect  in  situ  conditions  (Urmos  and 
Wilkens,  1993).  Indeed,  the  large  gradients  that  characterize  curves  1, 
2,  4,  and  5  are  consistent  with  in  situ  measurements  made  in  non- 

carbonate  muddy  sediments  (Fig.  19)  which  record  significant  velocity 

gradients  in  the  upper  20  m  of  sediment.  For  comparison,  curves  1  and 
4  (Fig.  18)  are  also  plotted  in  Figure  19. 

Compressional  Wave _ Atten.uaU.QP: 

The  most  comprehensive  compilations  of  attenuation  measure¬ 
ments  in  marine  sediments  have  been  assembled  by  Hamilton  (1972, 

1976a,  and  1987)  and  Kibblewhite  (1989).  The  acoustic  community 
continues,  for  the  most  part,  to  use  the  work  of  Hamilton  for  estimating 
attenuation  in  marine  sediments  (Kibblewhite,  1989).  Hamilton's  view, 
most  recently  stated  in  1987,  is  that  attenuation  (a)  varies  with  the 

first  power  of  frequency  (f)  and  can  be  expressed  as  a  function  of 
frequency  in  kilohertz: 


a(dB/m)  =  kP» 


where  k(dB/m/kHz)  is  a  constant  and  n  is  assumed  to  be  1.0.  Thus, 
compressional  wave  attenuation  can  be  calculated  for  any  particular 
frequency  using  the  above  equation. 

The  linear  frequency  dependence  of  attenuation  has  been 
questioned  (Stoll,  1977,  1985,  and  1989).  Moreover,  attenuation  versus 
frequency  relationships  indicate  basic  differences  between  the  response 
of  fluid-saturated,  coarse-and  fine-grained  sediments  (Stoll,  1985).  The 
sediments  within  the  Natural  Laboratory  area  are  assumed  to  be 
predominantly  fine-grained.  Therefore,  only  attenuation  data  relating 
to  fine-grained  sediments  were  taken  from  the  compilations  of  Hamilton 
(1972,  1976a,  and  1987)  and  used  to  generate  Figure  20.  Kibblewhite's 
review  (1989)  did  not  produce  any  additional  attenuation  data  relating 
to  fine-grained  sediments  and  a  search  of  1989  to  recent  literature 
produced  only  three  relevant  publications  (Chotiros  et  al.,  1989; 
Richardson  et  al.,  1992;  and  Rogers  et  al.,  1993).  Thus,  virtually  all 
published  compressional  wave  attenuation  measurements  (Table  4), 
relating  to  fine-grained  sediments,  appear  in  Figure  20. 

Of  particular  interest  in  Figure  20  is  the  nature  of  the  relationship 
between  attenuation  and  frequency.  A  curve  fit  of  the  data  shows  that 
attenuation  is  dependent  on  frequency  to  the  power  of  1.1115.  Over 
the  frequency  range  shown  the  attenuation  can  be  calculated  using  the 
following  equation: 


a(dB/m)  =  0.000025312  fl-lU5 

The  stippled  area  in  Figure  20  corresponds  to  the  range  of  predicted 
attenuation  values  for  "silts"  calculated  from  the  Biot  model  (Stoll,  1985 
and  1989).  It  can  be  seen  that  the  power  curve  fit  is  more  linear  in 
appearance  than  the  "silt"  envelope.  Nevertheless,  the  above  equation 
will  produce  values  that  are  generally  consistent  with  the  predicted 
range  of  attenuation  values  up  to  frequencies  of  about  100  kHz. 

Strickly  speaking,  the  data  in  Figure  20  do  not  support  the  linear 
dependence  of  attenuation  and  frequency  as  advocated  by  Hamilton 


Table  4.  Compratalonal  Wave  Attenuation  Data 


Rttarancf 

Publication 

Matarlal/Araa 

Oapth 

Fraquancy 

Attsnuatlon 

K 

Yaar 

(m) 

(HZ) 

(dB/m) 

(dB/m/kHz) 

Shumway 

mo 

day  ay  tit  (ava  12) 

0 

20500 

1  09 

0  053 

Shumway 

mo 

day  ay  sA  (ava  16) 

0 

22600 

1  67 

0  0732 

Shumway 

I960 

ctayay  sat  (ava. 20) 

0 

30400 

3.2 

0  105 

Wood  A  WMtOO 

1964 

mud 

0 

4000 

0  279 

0  0698 

Wood  A  Wsaton 

1964 

mud 

0 

6000 

0  479 

0  0595 

Wood  A  Waston 

1964 

mud 

0 

16000 

0  951 

0  0594 

Wood  A  Waston 

1964 

mud 

0 

32000 

1  8 

0  0564 

Wood  A  Wsston 

1964 

mud 

0 

46000 

3  02 

0  0663 

Bannstt 

I960 

abysr.4  plains  (ava) 

0 

12000 

0  7 

0  056 

Ulonska 

1968 

clay  ay  sik 

0 

250000 

26  2 

0  105 

McLaroy  A  DaLoach 

1968 

silty  ctay.(s«a  6) 

0 

15000 

1  08 

0  072 

McLaroy  A  DaLoach 

1968 

clay-stft.(srta  9) 

0 

15000 

1  12 

0  075 

McCann  A  McCann 

1969 

day 

0 

366000 

31  9 

0  087 

Barzon  #<  a< 

1969 

day 

140 

25 

0  0035 

0  139 

Schirroar 

1971 

mud  (day-silt) 

0 

1200 

0  29 

0  242 

Andaraon  A  Blackman 

1971 

calcaraous  day 

250 

200 

0  0031 

0  016 

Andaraon  A  Blackman 

1971 

caicaraous  day 

250 

315 

0  0128 

0  04 

Andaraon  A  Blackman 

1971 

calcaraous  day 

250 

250 

0  0069 

0  026 

Naprochnov 

1971 

taysn,  Bangs!  Bay 

350 

110 

0  0165 

0  15 

Naprochnov 

1971 

layan.  Arabian  Saa 

400 

110 

0.0156 

0  142 

Naprochnov 

1971 

layart.  S  indian  Ocn 

200 

110 

0  02 

0  182 

Naprochnov 

1971 

layan,  N.lndian  Ocn 

400 

110 

0  0182 

0  166 

Naprochnov 

1971 

layan,  Japan  Saa 

250 

110 

0  02 

0  182 

Naprochnov 

1971 

layan.  Black  Saa 

350 

110 

0.0185 

0  15 

Naprochnov 

1971 

Ob  T  ranch 

350 

110 

0  0174 

0  158 

Hamilton 

1972 

day  ay  si* 

0 

7000 

12 

0  171 

Hamilton 

1972 

day  ay  sit 

0 

14000 

2  4 

0  171 

Hamilton 

1972 

day  ay  sik 

0 

14000 

2  3 

0  164 

Hamilton 

1972 

ctayay  sik 

0 

14000 

1 

0  071 

Hamilton 

1972 

ctayay  sik 

0 

14000 

0  8 

0  043 

Hamilton 

1972 

ctayay  sik 

0 

23200 

2.8 

0  121 

Hamikon 

1972 

day  ay  aik 

0 

25000 

4 

0.16 

Hamilton 

1972 

ctayay  aik 

0 

1 00000 

16 

0  16 

Tyca 

1976 

carbons  ta 

28 

4000 

0.12 

003 

Dicua 

1976 

turbidita 

106 

65 

0  002 

0031 

Halmbargaf  at  al 

1979 

7  (Baring  Saa) 

150 

10 

0  000295 

0  0295 

Faria  at  al. 

1980 

mud 

0 

100 

0  0041 

0  041 

Faria  at  al. 

1960 

mud 

1 

250 

0017 

0  086 

Wroistad 

1980 

turbidita 

228 

80 

0  004 

0  05 

Chapman 

1960 

turbidita 

200 

160 

0  002 

0  013 

Mltchak  A  Focka 

1960 

turbidita 

0 

213 

00065 

0  04 

Mttchal  A  Focka 

1980 

turbidita 

0 

213 

0  0055 

0  026 

Rubano 

1980 

clay/silt 

0 

20 

0  0006 

0  03 

Rubano 

1980 

clay/silt 

0 

250 

0  0094 

0  036 

Hoopar  atai 

1981 

turbidita 

0 

163 

0.0054 

0  086 

Frisk  at  al 

1961 

slt/day 

80 

220 

0  0015 

0.007 

Tyca 

1961 

day  ay  sift? 

14 

4000 

0  29 

0  072 

Tyca 

1981 

ctayay  skt? 

16 

4000 

0.25 

0  062 

Tyca 

1981 

ctayay  a*!? 

18 

4000 

0  29 

0072 

Tyca 

1961 

hampaiagic  turbid 

27 

4 000 

0  26 

0.065 

Tyca 

1961 

hamipatagic  turbid. 

21 

4000 

0  19 

0  046 

Tyca 

1981 

hampatagic  turbid 

25 

4000 

0  28 

0.07 

Tyca 

1961 

hamipatagic  day 

14 

4000 

0  26 

0.065 

Tyca 

1961 

day  ay  silt? 

15 

4000 

0.21 

0  052 

Dicus  A  Andaraon 

1962 

turbidita 

0 

65 

0  0026 

0  04 

Dicus  A  Andaraon 

1982 

turbidita 

0 

65 

0  0033 

0  051 

Dicus  A  Andaraon 

1982 

caicaraous  day 

0 

65 

0  0176 

0.271 

Dicus  A  Andaraon 

1962 

turbidita 

0 

65 

0.0013 

0  02 

Tindla 

1962 

sifi/clay 

SO 

80 

0  0045 

0  056 

Tindla 

1982 

silt/clay 

50 

140 

0  0078 

0.056 

Stoll  A  Houtz 

1983 

sand/siivclay 

350 

60 

0  00756 

0  126 

Stoll  A  Houtz 

1983 

sand/sHt/clay 

350 

10 

0  00013 

0.013 

Stoll  A  Houtz 

1963 

sand/silvelay 

350 

20 

0  00113 

0.056 

Focka 

1983 

silt/clay 

0 

200 

0  006 

0  027 

Jacobson  at  al 

1984 

turbidita 

0 

95 

0  004 

0  042 

Baikar  A  Haimbargsr 

1986 

?  (Baring  Saa) 

150 

32 

0  0014 

0  045 

Barkar  A  Haimbargsr 

1986 

?  (Baring  Saa) 

150 

18 

0  0005 

0  027 

Frisk  at  al 

1986 

silt/clay 

26 

220 

0  0039 

0  018 

Frisk  at  al. 

1986 

silt/clay 

10 

220 

0  007 

0  032 

Duckworth  A  Baggaroor 

1966 

turbidita 

100 

27  5 

0  0019 

0  069 

Chotiros 

1989 

clay  (ava.26) 

30 

400000 

39  1 

0  098 

Richardson  at  al. 

1992 

Hotocanc  days 

0 

400000 

59 

0.15 

Richardson  at  al 

1992 

Hampaigic  sky  clays 

0 

400000 

50 

0  12 

Rogara  at  al. 

1993 

sand  A  siky  day 

100 

50 

0  000557 

0.01 1 1 

Rogaraatal. 

1993 

sand  A  siky  day 

100 

75 

0  000753 

0.01 

Rogara  at  al. 

1993 

sand  A  siky  day 

100 

175 

0 00213 

0  0123 

Rogaraatal 

1993 

sand  A  siky  day 

100 

275 

0  00381 

0  0131 

Rogara  at  al 

1993 

sand  A  siky  clay 

100 

375 

0  00524 

0  014 

Rogaraatal. 

1993 

sand  A  siky  day 

100 

525 

0  00802 

0  0153 

Rogara  at  al. 

1993 

sand  A  siky  day 

100 

600 

0 00955 

0  0159 
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(1987),  even  for  fine-grained  sediments.  Nevertheless,  the  relationship 
is  close  enough  to  being  linear  that  a  first  order  approximation  of 
attenuation  may  be  acceptable  for  purposes  of  acoustic  modeling.  A 
linear  fit  to  the  data  produces  that  following  equation: 

a(dB/m)  =  -0.027116  +  0.00010532x 

For  most  of  the  frequency  range  shown  in  Figure  20  this  equation 
provides  attenuation  values  that  are  very  close  to  the  values  derived 
from  the  power  law  equation,  but  it's  use  is  not  recommended.  At  500 
Hz  both  equations  give  essentially  the  same  attenuation  value  (i.e.,  the 
curves  intersect).  Below  500  Hz,  however,  the  values  derived  from  the 
linear  equation  decrease  rapidly  (i.e.,  the  linear  curve  drops  rapidly 
away  from  the  power  curve);  below  260  Hz  they  become  negative.  At 
frequencies  above  100  kHz  the  power  equation  gives  attenuation  values 
that  increase  at  a  much  faster  rate  than  values  derived  from  the  linear 
equation  (i.e.,  the  curves  diverge  rapidly). 

In  Figure  21  the  same  attenuation  data  used  for  Figure  20  are 
expressed  as  the  attenuation  coefficient  k  and  plotted  against  depth. 
Superimposed  on  the  data  are  attenuation-depth  profiles  from  Hamilton 
(1976a),  Mitchell  and  Focke  (1980),  Jacobson  et  al.  (1981),  and  Brienzo 
(1992).  For  the  present,  these  profiles  appear  to  be  the  only  choices 
that  modelers  have  for  predicting  compressional  wave  attenuation  at 
significant  depths  in  the  sediment  column.  Although  Hamilton’s  profile 
is  widely  used,  it  can  seen  that  his  silt-clay  profile  is  based  on  sparse, 
scattered  data  (Fig.  22).  When  plotted  with  additional  data,  as  in  Figure 
21,  the  curve  fit  is  not  improved.  Jacobson's  profile  is  based  on 
attenuation  values  for  a  region  of  thick,  sand/silt  turbidite  deposits  (i.e., 
Bengal  deep-sea  fan)  and,  as  such,  may  not  be  suitable  for  modeling 
most  deep,  open-ocean  areas.  Brienzo's  profile  is  also  from  an  area  of 
turbidite  deposition  (Monterey  Fan),  yet  the  attenuation  values  are 
significantly  lower  than  Jacobson's.  Indeed,  they  are  very  close  to 
values  represented  by  the  Mitchell-Focke  profile  that  was  derived  from 
bottom  loss  data  taken  in  a  region  of  calcareous  ooze/calcareous  clay 
sediment.  Considering  sediment  type,  the  Mitchell-Focke  profile 
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represents  a  more  appropriate  choice  for  geoacoustic  modeling  in  deep 
ocean  areas.  The  following  regression  equation, 

K  =  0.0262  +  8.301xl0*5d  +(-2.483xl0-7)d2  +  1.455xl0*7d3 

where  d  =  depth  and  K  is  in  dB/m/kHz,  will  essentially  reproduce  their 
profile  (for  the  upper  500  m  of  sediment).  Matthews  (1982)  notes  that 
the  Mitchell-Focke  attenuation  values  yield  plane  wave  bottom  loss 
calculations  in  close  agreement  with  measured  bottom  loss  values. 

Shear _ Wave _ Attenuation: 

Very  little  information  is  available  regarding  the  variations  of 
shear  wave  attenuation  with  depth,  particularly  in  unlithified,  high- 
porosity,  fine-grained  sediments  (Table  5).  For  modeling  purposes 
Hamilton  (1976c,  1987)  has  assumed  that  shear  wave  attenuation 

varies  with  depth  in  the  sea  floor  proportionally  with  compressional 
wave  attenuation  according  to  the  following  equation: 

KSz*(KSo/kpo)  Kpz 

Accordingly,  shear  wave  attenuation  can  be  expected  to  track  the 
curves  shown  in  Figure  21  but  at  much  higher  K  values  (Hamilton 
suggests  a  surface  value  of  K  =  17.3  for  unconsolidated  clays).  This 
assumption,  however,  does  not  hold  with  more  recent  data. 

Sauter  (1987)  used  spectral  ratios  of  Scholte  waves  to  estimate 
1/Q  in  the  sediments  for  a  deep-water  site  off  the  California  coast.  His 
work  represents  a  well-constrained  study  detailing  changes  in  shear 
wave  attenuation  within  the  upper  60  m  of  sediment.  In  Figure  23 
Sauter’s  values  for  1/Q  have  been  converted  to  K(dB/m/kHz)  and 
plotted  with  depth  along  with  various  shear  wave  attenuations  reported 
by  other  investigators.  Sauter’s  study  indicates  a  much  higher  shear 
wave  attenuation  at  the  sediment  surface  (k  =  40)  than  Hamilton 
(1976c)  suggests,  and  also  shows  that  the  attenuation  decreases  rapidly 
within  the  upper  5  m  of  sediment.  Between  5-20  m  the  change  is  more 
gradual,  and  from  20-60  m  the  attenuation  is  practically  constant. 
Shear  wave  estimates  reported  by  Jensen  and  Schmidt  (1986),  Carter  et 
al.  (1986),  and  Bromirski  et  al.  (1992)  for  depths  of  60  and  180  m 
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Tabla  S.  Shaar  Wava  Attanuation  Data 


Ratoranca 

Publication 

Malar  lal/Araa 

Dapth 

Frequency 

Attanuation 

K 

Year 

("»> 

(Hi) 

(dB/m) 

(dB/m/kHz) 

Bucket  at  al. 

1964 

sand 

0 

20 

0.328 

16.4 

Bucker  at  al. 

1964 

sand 

0 

25 

0.197 

7.88 

Wamck 

1974 

soft  mud 

5 

5 

0.095 

18.9 

Joynar  at  al. 

1976 

aofl  mud 

0 

5 

0.284 

56  8 

Schirmer 

1980 

North  Saa  (130m) 

15 

4.5 

0.007 

1.56 

Holt,  at  al. 

1983 

sand 

2 

35 

0.6 

17.1 

Brocher  at  al. 

1983 

sand 

7 

5 

0.0038 

0.769 

Janaan  A  Schmidt 

1986 

sand/silt 

0 

3 

0.016 

5  29 

Jensen  &  Schmidt 

1986 

sand/silt 

60 

3 

0.002 

0.5 

Rauch 

1986 

soli  sadimant  (mud?) 

15 

3.5 

0.0056 

1.6 

Cartar  at  al. 

1986 

biosillcaous  clay 

180 

18 

0.0205 

1.137 

Sautar 

1987 

clay/ooza 

12 

4 

0.01 

2.43 

Sautar 

1987 

clay/ooza 

14 

4 

0.0085 

2.13 

Sautar 

1987 

clay/ooza 

10 

4 

0.011 

2.73 

Sautar 

1987 

clay/ooza 

6 

4 

0.02 

5.05 

Sautar 

1987 

clay/ooza 

8 

4 

0.014 

3.59 

Sautar 

1987 

clay/ooza 

40 

4 

0.0033 

0.82 

Sautar 

1987 

clay/ooza 

60 

4 

0.0023 

0.57 

Sautar 

1987 

clay/ooza 

20 

4 

0.0064 

1.61 

Sautar 

1987 

clay/ooza 

16 

4 

0.0076 

1.9 

Sautar 

1987 

clay/ooza 

t8 

4 

0.0071 

1.76 

Sautar 

1987 

clay/ooza 

1 

4 

0.112 

28.06 

Sautar 

1987 

clay/ooza 

0.5 

4 

0.136 

34.11 

Sautar 

1987 

clay/ooza 

0 

4 

0.164 

40  97 

Sautar 

1987 

clay/ooza 

2 

4 

0.091 

22.76 

Sautar 

1987 

clay/ooza 

3.5 

4 

0.039 

9.74 

Sautar 

1987 

clay/ooza 

3 

4 

0.055 

13.64 

Sautar 

1987 

clay/ooza 

2.5 

4 

0.072 

17.95 

Sarano 

1990 

zaolltic  Clay 

35 

15 

0.118 

7.84 

Muir  at  al. 

1991 

mud 

2.5 

4 

0.017 

4.22 

Muir  at  al. 

1991 

mud 

7.5 

4 

0.008 

2.08 

Mur  at  al. 

1991 

mud 

5 

4 

0.012 

3.01 

Muir  at  al. 

1991 

mud 

30 

4 

0.0047 

1.17 

Muir  at  al. 

1991 

mud 

15 

4 

0.0054 

1.34 

Muir  at  al. 

1991 

mud 

0 

4 

0.0345 

6.64 

Bromirski  at  al. 

1992 

biosiliceous  clay 

180 

18 

0.0025 

1.4 

Bromirski  at  al. 

1992 

biosillcaous  Clay 

180 

3 

0.0042 

1.4 

15 


(Table  5),  coincide  with  Sauter’s  profile  and  support  the  observation 
that  attenuation  is  essentially  constant  below  20  m. 

Decreasing  shear  wave  attenuation  in  the  upper  4  m  of  muddy 
sediment  (Fig.  24)  is  also  documented  by  Muir  et  al.  (1991)  from 
measurements  made  in  short  (6  m)  boreholes.  Interestingly,  their 
results  indicate  that  attenuation  varies  with  frequency  to  the  power  of 
1.5,  contradicting  Hamilton’s  (1987)  recommendation  that  shear  wave 
attenuation  be  considered  to  vary  linearly  with  frequency.  It  is 
apparent  in  Figure  24  that  the  attenuation  curves  for  the  different 
frequencies  converge  in  the  upper  4  m  of  sediment.  Figure  25  combines 
the  data  in  Figure  24  with  the  data  in  Figure  23  for  the  upper  20  m  of 
sediment  and  shows  that  shear  attenuation  at  low  frequencies  (100  Hz 
and  less)  is  most  variable  in  the  first  4-5  m  of  sediment  and  is 
essentially  constant  below  about  12  m.  A  comparison  (Fig.  26)  of 
Sauter’s  data  (1987)  with  Scholte  wave  data  presented  by  Muir  et  al. 
(1991)  also  shows  large  differences  in  attenuation  in  the  upper  5  m  of 
sediment  (due  principally  to  different  sediment  types  -  see  Table  5)  and 
convergence  of  the  data  below  about  20  m. 

RESULTS: 

The  geoacoustic  parameters  (porosity,  wet  bulk  density, 
compressional  wave  velocity  and  attenuation,  shear  wave  velocity  and 
attenuation)  discussed  in  the  previous  section  are  presented  below  in 
the  form  of  a  geoacoustic  model  that  incorporates  the  variations  of  each 
parameter  to  a  depth  of  300  m  (estimated  sediment  thickness  for  a 
large,  midocean  ridge  flank  sediment  pond). 


GEOACOUSTIC  MODEL: 


Material  Depth  Vp  Vs  p  n  KpKs 

Type  (mbsf)  (m/sec)  (m/sec)  (g/cm3)  (%)  (dB/m-kHz) 


Nannofossil 

clay/ooze 


Basalt 


0 

1487 

23 

1.574 

68.1 

0.0266 

41.0 

1 

1488 

67 

1.575 

68.1 

0.0267 

28.1 

2 

1490 

80 

1.576 

68.0 

0.0267 

22.8 

5 

1494 

110 

1.579 

67.9 

0.0268 

6.93 

10 

1502 

145 

1.583 

67.7 

0.0270 

2.73 

20 

1516 

180 

1.592 

67.3 

0.0275 

1.61 

30 

1531 

208 

1.600 

66.9 

0.0278 

1.21 

40 

1545 

230 

1.610 

66.4 

0.0282 

0.82 

50 

1559 

249 

1.618 

66.0 

0.0285 

0.69 

60 

1572 

266 

1.627 

65.6 

0.0292 

0.57 

80 

1599 

294 

1.644 

64.8 

0.0301 

0.57 

100 

1624 

319 

1.662 

63.9 

0.0310 

0.57 

120 

1649 

340 

1.679 

63.1 

0.0315 

0.57 

140 

1672 

359 

1.697 

62.2 

0.0322 

0.57 

160 

1695 

376 

1.714 

61.4 

0.0327 

0.57 

180 

1717 

392 

1.732 

60.6 

0.0331 

0.57 

200 

1738 

407 

1.749 

59.8 

0.0335 

0.57 

220 

1759 

421 

1.767 

58.9 

0.0337 

0.57 

240 

1778 

434 

1.784 

58.1 

0.0340 

0.57 

260 

1797 

447 

1.802 

57.2 

0.0342 

0.57 

280 

1816 

459 

1.819 

56.4 

0.0342 

0.57 

300 

1833 

470 

1.837 

55.6 

0.0338 

0.57 

4100 

2100 

2.810 

0.0200 

0.07 

Notes: 

1.  Depth: 

mbsf  =  meters  below  seafloor 
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2.  Vp  ('Compressional  wave  velocity): 


Vp  =  compressional  wave  velocity  =  1.487  +  0.00013369d  +  AV  p 
where  d  =  depth(m)  and  AVp  =  0.66  (l-e-°  00208d). 

The  first  part  of  the  above  equation  is  the  same  as  the  equation 
shown  in  Figure  15  except  that  the  initial  velocity  at  d  =  0  m  has 
been  changed  from  1.529  km/sec  to  1.487  km/sec.  In  situ 
measurements,  even  in  clayey  silts,  record  surficial  velocities  below 
1500  m/sec  (Hamilton,  1972).  The  initial  velocity  of  1.487  km/sec 
is  taken  from  an  equation  derived  by  Lovell  and  Ogden  (1984)  for 
surficial  carbonate  samples. 

AVp  is  a  correction  factor  derived  by  Urmos  and  Wilkens  (1993)  to 
compensate  for  differences  in  log  (in  situ)  and  laboratory  velocities 
measured  in  calcareous  sediments  having  carbonate  contents  of  60% 
or  greater.  The  majority  of  the  samples  compiled  for  this  report 
have  carbonate  contents  greater  than  60%  (see  Figure  1  and  Table 

3). 

3.  Vs  fShear  wave  Velocity): 

Vs  =  shear  wave  velocity  =  62.5^0-3537  (from  Lovell  and  Ogden, 
1984) 

Vs  =  23.0  m/sec  at  d  =  0  m  is  an  average  value  derived  from 
measurements  reported  by  Richardson  (1987;  Table  1)  for  fine¬ 
grained  (noncalcareous)  sediments.  Vs  =  67  m/sec  at  d  =  1  m  is 
taken  from  the  (overlapping)  Muir,  Stoll,  Lewis,  and  Davis  curves  in 
Figure  19.  Because  these  four  curves  coincide  closely  with  curve  4 
in  Figure  18  (See  Fig.  19),  the  equation  representing  curve  4  (above) 
is  used  to  compute  Vs  for  2-300  m. 

4.  p  (Bulk  Density): 


p  =  1.5743  +  0.0008742d,  where  d  =  depth(m).  This  equation  is 
derived  from  a  linear  curve  fit  of  the  data  in  Figure  16. 


5.  n  (Porosity): 

n  =  143.46  +  (-47.85 lp),  where  p  =  bulk  density.  This  equation  is 
derived  from  a  linear  curve  fit  of  the  data  in  Figure  4.  The  equation 
shown  in  Figure  17  will  produce  comparable  values. 

6.  Kp  (compressional  wave  attenuation): 

Kp  =  0.0266  +  6.9  xl0*5d  +  (-1.5  xl0'7)d2,  where  d  =  depth(m) 

The  above  equation  represents  the  Mitchell  and  Focke  curve  shown 
in  Figure  21.  The  equation  does  not  hold  below  500  m. 

7.  Ks  (shear  wave  attenuation): 

Ks  values  are  picked  directly  from  the  attenuation  model  generated 
by  Sauter  (1987)  and  shown  in  Figures  26  and  27.  Sauter's  model  is 
derived  from  measurements  made  at  a  deep  water  site  (3800  m) 
near  DSDP  Site  469.  Although  the  upper  42  m  of  sediment  in  this 
area  is  described  as  a  clay  with  a  carbonate  content  less  than  30% 
(Yeats,  et  al.,  1981),  it  is  thought  to  be  the  best  choice,  at  this  time, 
for  estimating  attenuation  as  a  function  of  depth  in  soft  sediments 
regardless  of  carbonate  content. 


MODEL  APPLICABILITY: 

The  model  presented  in  this  report  is  generic  and  should  be 
applicable  to  most  areas  that  are  covered  with  soft,  fine-grained, 
calcareous  sediment.  It  is  important  to  note,  however,  that  high- 
resolution,  deep-towed  multichannel  seismic  data  (Gettrust  et  al.,  1988; 
Rowe  and  Gettrust,  1989;  Rowe  and  Gettrust,  1992;  Rowe  and  Gettrust, 
1993)  resolve  significant  lateral  and  vertical  variability  in 
compressional  velocity  over  distances  of  about  ten  to  one  hundred 
meters.  Fine-scale  spatial  variability  such  as  this  is  not  addressed  by 
presently  accepted  geoacoustic  modeling  methods  (Hamilton,  1972, 
1976a,  1976b,  1987;  Hamilton  and  Bachman  1982)  nor  by  the  model 
presented  here.  Moreover,  these  data  reveal  that  compressional 
velocity  gradients  in  the  upper  few  hundred  meters  of  the  sediment 
column  can  be  substantially  lower  than  conventional  geoacoustic  models 
predict. 
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